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ABSTRACT 

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova 
survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CL012Car ( z = 1.28), 
SN CLN12Did (z = 0.85), and SN CLAllTib (z — 1.14), are located behind three different clusters, 
MACS J1720. 2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and A383 (z = 0.187), respectively. Each 
SN was detected in Hubble Space Telescope optical and infrared images. Based on photometric classification, we 
find that SNe CL012Car and CLN12Did are likely to be Type la supernovae (SNe la), while the classification 
of SN CLAllTib is inconclusive. Using multi-color light-curve fits to determine a standardized SN la luminosity 
distance, we infer that SN CL012Car was ~1.0 =b 0.2 mag brighter than field SNe la at a similar redshift and 
ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is ~0.2 =b 0.2 mag brighter than field 
SNe la. We derive independent estimates of the predicted magnification from CLASH strong+weak-lensing maps 
of the clusters (in magnitude units, 2.5 log 10 /x): 0.83 =b 0.16 mag for SN CL012Car, 0.28 =b 0.08 mag for SN 
CLN12Did, and 0.43 d= 0.11 mag for SN CLAllTib. The two SNe la provide a new test of the cluster lens model 
predictions: we find that the magnifications based on the SN la brightness and those predicted by the lens maps 
are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe la (from 
the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct 
determination of absolute magnifications. 

Key words: cosmology: observations - galaxies: clusters: general - gravitational lensing: weak - 
supernovae: general 
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1. INTRODUCTION 

Supernovae (SNe) have proven to be important tools for 
studying the universe. Type la SNe (SNe la) have been in- 
strumental in the discovery that the expansion of the universe is 
accelerating, and probing the dark energy driving the accelera- 
tion (e.g., Riess et al. 1998; Perlmutter et al. 1999; Wood-Vasey 
et al. 2007; Kessler et al. 2009a; Hicken et al. 2009; Sullivan 
et al. 2011; Suzuki et al. 2012 and references therein). Core- 
collapse (CC) SNe trace the star-formation history of the uni- 
verse (e.g., Dahlen et al. 2004; Bazin et al. 2009; Botticella et al. 
2012). Gravitational lensing of SNe can augment their cosmo- 
logical utility, in both the strong- and weak-lensing regimes. 
The seminal paper of Refsdal (1964), which considered mea- 
suring galaxy masses and the Hubble parameter from lensed 
image time delays, in fact envisaged SNe as the background 
sources. 35 While no multiply imaged SN has yet been discov- 
ered, several applications of lensed SNe are currently being 
explored. 

Gravitational lensing allows massive galaxy clusters to be 
used as cosmic telescopes, increasing the apparent brightness 
of distant sources that would otherwise be too faint to detect 
(Sullivan et al. 2000; Kneib et al. 2004; Stanishev et al. 2009). 
Ground-based near-infrared (IR) surveys of galaxy clusters are 
being used to find lensed SNe, including a CC SN with redshift 
z — 1.7 (Goobar et al. 2009; Amanullah et al. 2011 and 
references therein). Lensed SNe la are particularly useful, as 
their standardizable luminosities allow for a direct measurement 
of the absolute lensing magnification (Kolatt & Bartelmann 
1998). 

The Cluster Lensing And Supernova survey with Hubble 
(CLASH; Postman et al. 2012) has observed 25 galaxy clusters 
in 16 Hubble Space Telescope (. HST) broadband filters covering 
the near-ultraviolet (UV) to the near-IR with the Advanced 
Camera for Surveys (ACS) and the Wide Field Camera 3 
(WFC3). The WFC3 infrared channel (WFC3/IR) is especially 
important for finding and studying high-redshift SNe, efficiently 
surveying large sky areas for SNe whose peak flux has redshifted 
from the optical to the IR. The SN discoveries and follow-up 
observations from CLASH are discussed at length by Graur 
et al. (2014). Here we focus on SNe behind CLASH clusters (in 
the “prime” fields centered on the clusters); these have not been 
included in the Graur et al. (2014) “parallel” field sample. 

All distant SNe are lensed (magnified, or more often, demag- 
nified) to some extent by large-scale structure along the line 
of sight. This adds a statistical “nuisance” in cosmological in- 
ferences from SNe la that can be corrected either assuming an 
overall magnification distribution (e.g., Holz & Linder 2005; 
Martel & Premadi 2008) or with convergence estimates along 
each SN light of sight using simplified scaling relations applied 
to the nearby foreground galaxies observed (Kronborg et al. 
2010; Jonsson et al. 2010; Smith et al. 2014, and references 
therein). 

Individual SNe that are detectably magnified have been found 
only rarely in recent wide-field surveys. A gravitational lens has 
to be precisely aligned with the background source to produce 
large magnifications, and this is significantly rarer for SNe 


34 Hubble Fellow. 

35 Refsdal (1964) also presciently concluded with the idea that in addition to 
lensed SNe, “star-like objects with intense emission both in the radio range 
and optical range” that “have been recently discovered” may also be valuable 
lensed sources, “giving a possibility of testing the different (cosmological) 
models.” 


compared to background galaxies as sources. Their short-lived 
durations also make lensed SNe rarer than even lensed quasars 
(Oguri & Marshall 2010). Moreover, the survey footprints of the 
Supernova Legacy Survey (SNLS) and the SDSS-II Supernova 
Survey did not contain many massive galaxy clusters to act as 
lenses for background SNe (Graham et al. 2008; Dilday et al. 
2010 ). 

Kronborg et al. (2010) inferred that the most magnified SN la 
in SNLS had a magnification factor (i = 1.27, corresponding to 
Ara^ = 2.5 log 10 /x = 0.26 mag, calculated from the location 
and photometry of foreground galaxies along the line of sight. 
Jonsson et al. (2010) came to a similar conclusion, finding that 
the SNLS sample did not probe SNe la with A 0.25 mag 
(however, see Karpenka et al. 2013, who find little to no lensing 
of the SNLS sample). Given the typical ~0.2 mag scatter in 
standardized SN la luminosities, a lensing signal at this level 
cannot be clearly ascertained for an individual object, although 
it can be detected statistically for a larger sample. 

Lensing magnification by foreground galaxies is expected to 
play a larger role for more distant SNe. The Hubble Deep Field 
SNIal997ffaU « 1.7 (Gilliland etal. 1999; Riess etal. 2001) is 
calculated to have a Am^ = 0.34 =b 0. 12 mag from lens galaxies 
projected nearby (Benitez et al. 2002). The uncertainties in the 
SN light curve and cosmological model preclude a definitive 
detection of the expected magnification signal. Rather, in this 
case the estimated lensing magnification has been applied as a 
correction to the inferred SN luminosity distance to constrain 
cosmological parameters, though care should be taken that this 
is done self-consistently (Zitrin et al. 2013). 

Supernova properties were used by Quimby et al. (2013) to 
posit that the recent SN PSl-lOafx at z = 1.39 was an SN la 
magnified by a factor of 30 (A = 3.72 ± 0.18 mag)! 
However, there is no lens observed in this system, and the 
authors conclude that it must be dark. Chornock et al. (2013) 
offer an alternate model, that this SN is best explained as a 
superluminous CC object, without significant magnification. 

Ideally, we would like to confirm the detection of gravitational 
lensing magnification of an SN by independently checking the 
magnifications derived from the SN brightness and a lens-model 
prediction for consistency. In this paper, we discuss three grav- 
itationally lensed SNe (SN CL012Car, SN CLN12Did, and SN 
CLAllTib) behind CLASH clusters. We present photometric 
observations of the SNe and spectroscopy of their host galax- 
ies in Section 2. In Section 3, we discuss classification of the 
SNe and model light-curve fits. The gravitational magnifications 
from lens models are derived in Section 4, and we confront 
the models against the observations and discuss the results in 
Section 5. The following cosmological parameters are adopted 
throughout the paper (Sullivan et al. 2011): Ho = 71.6 km s -1 
Mpc -1 , = 0.27, Qa = 0.73, and = 0. 

2. OBSERVATIONS 
2.7. SNCLOUCar 

SN CL012Car was discovered in the CLASH images of 
MACSJ1 720.2+3536 (Ebeling et al. 2010). The SN was detected 
in 2012 June with both ACS and WFC3/IR. Figure 1 shows 
the location of the SN (a = 17 h 20 m 21 s .03, 5 = +35 o 36 / 40?9, 
J2000), (L35 from the center of its host galaxy. The host is 
relatively bright in the IR (F160W = 21.50 zb 0.01 AB mag), 
but is significantly fainter in the optical images (F606W = 
24.21 =b 0.06 AB mag). 
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UT Date 

MJD 

Instrument 

Filter 

Exposure Time 

(s) 

Magnitude 

2012 May 5 

56052.58 

ACS 

F850LP 

991.0 

>25.546 

2012 May 5 

56052.64 

WFC3/IR 

F160W 

1408.8 

>25.347 

2012 May 9 

56056.02 

WFC3/IR 

F105W 

1408.8 

>26.016 

2012 May 9 

56056.04 

WFC3/IR 

F140W 

1005.9 

>25.808 

2012 May 22 

56069.68 

ACS 

F814W 

1007.0 

>26.666 

2012 Jun 17 

56095.66 

ACS 

F850LP 

1032.0 

24.597 ±0.130 

2012 Jun 17 

56095.68 

ACS 

F814W 

975.0 

25.148 ±0.082 

2012 Jun 17 

56095.73 

WFC3/IR 

F160W 

1408.8 

23.451 ±0.050 

2012 Jun 17 

56095.75 

WFC3/IR 

FI 10W 

1005.9 

23.856 ± 0.024 

2012 Jul 2 

56110.09 

WFC3/IR 

F105W 

1005.9 

23.750 ± 0.042 

2012 Jul 2 

56110.15 

WFC3/IR 

F140W 

1005.9 

23.271 ±0.025 

2012 Jul 16 

56124.09 

WFC3/IR 

F105W 

1005.9 

23.791 ±0.054 

2012 Jul 16 

56124.15 

WFC3/IR 

F140W 

1005.9 

23.250 ± 0.025 

2012 Jul 23 

56131.31 

WFC3/IR 

F160W 

455.9 

23.340 ± 0.063 

2012 Jul 23 

56131.38 

WFC3/IR 

F105W 

455.9 

24.138 ±0.056 


Notes. All supernova photometry is reported in Vega magnitudes. The nondetections are reported as 3cr upper limits. 



Figure 1. Top: CLASH HST composite false-color image of SN CL012Car. 
In this image, the red channel is an average of the WFC3/IR F105W, F110W, 
F125W, F140W, and F160W data. The green channel combines ACS F606W, 
F625W, F775W, F8 14W, and F850LP observations, and the blue channel is from 
ACS F435W and F475W data. The location of the SN is marked with a white 
circle. Bottom: the left panel shows a 5" square cutout of the WFC3/IR F160W 
image from 2012 June 17, with SN CL012Car marked by the arrow; the right 
panel shows the difference image after template subtraction. 

(A color version of this figure is available in the online journal.) 

We present aperture photometry for SN CL012Car in 
Table 1. These data include publicly available follow-up ob- 
servations collected in 2012 July by HST program GO- 12360 


(PI: S. Perlmutter). The flux from SN CL012Car was measured 
in difference images, using templates constructed with CLASH 
observations from 2012 March and April, more than 30 days 
before the first detection. For the ACS images, we used a fixed 
4 pixel (O'/ 20) radius aperture, whereas for WFC3/IR, the pho- 
tometry was measured in a 3 pixel radius (0'/27) aperture. 

We used the DEep Imaging Multi-Object Spectrograph 
(DEIMOS; Faber et al. 2003) on the Keck-II 10 m Telescope to 
obtain a ~2 hr spectrum of the host galaxy of SN CL012Car 
(2 x 1800 s, 1 x 1700 s, and 1 x 1500 s), with ~4 A resolution, 
covering the wavelength range 5000-9000 A. The spectrum 
showed evidence of a weak emission line located between two 
night-sky lines near 8506 A. For confirmation, we obtained a 
2.5 hr (5 x 1800 s) Gemini North Multi-Object Spectrograph 
(GMOS; Hook et al. 2004) nod-and-shuffle spectrum of the host, 
spanning 7460-9570 A with a resolution of ~2. 1 A. The GMOS 
spectrum also shows evidence for the faint line (marginally re- 
solved as a doublet), which we identify as [O ii] AA3726, 3728 
emission at z = 1.281 ± 0.001. A composite of the DEIMOS 
and GMOS spectra is shown in the top panel of Figure 2. The 
spectroscopic redshift for the CL012Car host galaxy lies near 
the peak of its photometric redshift (photo-z) probability den- 
sity function (PDF), but the photo-z PDF is broad, only limiting 
the redshift to 0.9 z 1.9. A full description of the CLASH 
photo-z estimation methods is given by Jouvel et al. (2014) and 
A. Molino et al. (in preparation). At z = 1.281, a separation of 
(X/35 from SN to host-galaxy center corresponds to a projected 
physical separation of ~3 kpc (modulo lensing magnification). 

We also analyzed the follow-up HST WFC3/IR G102 and 
G141 grism (Dressel 2012) data for SN CL012Car and its host 
galaxy (taken while the SN was still visible; program GO- 12360, 
PI: S. Perlmutter). The G102 spectrum has a dispersion of 25 A 
pixel -1 , covering 8000-11500 A, while the G141 spectrum has 
a dispersion of 47 A pixel -1 , covering 11000-17000 A. The 
exposure time with each grating was about 1 .3 hr, but this yielded 
only a weak signal, and we were unable to extract a clear SN 
spectrum. In these slitless grism observations, spectral features 
are effectively convolved with the spatial profile of the source, 
making sharp emission lines broader and less pronounced. There 
is an emission line in the G141 grism spectrum at 14963 A, 
consistent with Ha A.6563 at z ^ 1.28, but this is only a marginal 
detection given the low resolution of the data. 
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Figure 2. Spectra of the host galaxies of the three lensed SNe. Each panel shows observed wavelength along the bottom axis and rest wavelength on top. For display 
purposes, all spectra have been smoothed with a running 3 pixel median. Top: the host galaxy of SN CL012Car, at z = 1.281 ± 0.001. The plot shows the mean 
spectrum combining observations with Gemini+GMOS and Keck+DEIMOS. The single emission line at 8500 A is identified as the [O n] AA3726, 3728 doublet 
based on the photo-z and asymmetry of the line profile. Middle: the CLN12Did host galaxy, combining spectra from Gemini+GMOS and LBT+MODS. We derive 
^ = 0.851 ± 0.002 from cross-correlation with galaxy templates and identification of Ca and H absorption features, as indicated. Bottom: the left panel shows the 
spectrum of the SN CLAllTib host galaxy, observed with VLT+VIMOS. We derive z = 1.143 ± 0.001 from multiple absorption features and [C hi] A 1909 emission. 
This confirms the [O ii] emission-line redshift from the SALT+RSS spectrum shown on the right. 


Table 2 

SN CLN12Did Photometry 


UT Date 

MJD 

Instrument 

Filter 

Exposure Time 

(s) 

Magnitude 

2012 Feb 3 

55960.64 

ACS 

F625W 

1032.0 

26.699 ±0.171 

2012 Feb 3 

55960.65 

ACS 

F850LP 

1017.0 

25.066 ±0.150 

2012 Feb 18 

55975.67 

ACS 

F775W 

1032.0 

23.041 ± 0.016 

2012 Feb 18 

55975.69 

ACS 

F606W 

998.0 

24.186 ±0.017 

2012 Mar 3 

55989.84 

ACS 

F814W 

1032.0 

23.175 ±0.015 

2012 Mar 3 

55989.86 

ACS 

F435W 

1018.0 

>28.560 

2012 Mar 3 

55990.24 

ACS 

F625W 

1032.0 

24.182 ±0.023 

2012 Mar 3 

55990.25 

ACS 

F850LP 

1017.0 

22.929 ± 0.026 

2012 Mar 3 

55990.31 

ACS 

F475W 

1032.0 

27.057 ±0.138 

2012 Mar 3 

55990.32 

ACS 

F775W 

1013.0 

23.034 ± 0.016 

2012 Mar 3 

55990.37 

WFC3/IR 

FI 10W 

1508.8 

22.915 ±0.013 

2012 Mar 3 

55990.39 

WFC3/IR 

F160W 

1005.9 

23.041 ± 0.041 

2012 Mar 16 

56002.61 

ACS 

F606W 

1032.0 

25.491 ±0.038 

2012 Mar 16 

56002.62 

ACS 

F814W 

984.0 

23.410 ±0.018 

2012 Mar 16 

56002.67 

WFC3/IR 

F105W 

1305.9 

23.138 ±0.019 

2012 Mar 16 

56002.69 

WFC3/IR 

F140W 

1305.9 

23.540 ± 0.034 

2012 Mar 18 

56004.74 

ACS 

F475W 

1032.0 

>28.924 

2012 Mar 18 

56004.75 

ACS 

F850LP 

1001.0 

23.359 ± 0.037 

2012 Mar 29 

56015.52 

ACS 

F435W 

1032.0 

>28.625 

2012 Mar 29 

56015.53 

ACS 

F814W 

1017.0 

24.148 ±0.029 

2012 Apr 12 

56029.62 

ACS 

F850LP 

1032.0 

24.448 ± 0.079 

2012 Apr 12 

56029.63 

ACS 

F814W 

985.0 

25.072 ± 0.058 

2012 Apr 12 

56029.68 

WFC3/IR 

F160W 

1508.8 

23.316 ±0.038 

2012 Apr 12 

56029.70 

WFC3/IR 

FI 10W 

1005.9 

23.853 ± 0.038 


Notes. All SN photometry is reported in Vega magnitudes. The nondetections are listed as 3cr upper limits. 


2.2. SNCLN12Did a = 15 h 32 m 59 s .25, S = +30°21'42"8, J2000), and we present 

aperture photometry of the SN in Table 2. The SN was present in 
the first epoch of observations in each of the ACS and WFC3 /IR 
filters, so we did not have an SN-free template for subtraction. 
Nonetheless, the SN is separated by from its putative host 


SN CLN12Did was discovered in early 2012 in CLASH 
ACS and WFC3/IR images of RXJ1532.9+3021 (Ebeling et al. 
1998). Figure 3 shows the position of the SN (located at 
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Figure 3. Top: HST composite false-color image of SN CLN12Did and 
surroundings. The red channel is an average of WFC3/IR F105W, F110W, 
F140W, and F160W data from the CLASH program, the green channel combines 
the ACS F606W, F625W, F775W, F814W, and F850LP observations, and the 
blue channel is from ACS F435W and F475W data. The SN position is marked 
with a white circle, in the outskirts of its early-type host galaxy to the southeast. 
Bottom: the left panel shows a 5" square cutout of the ACS F850LP image from 
2012 March 3, with SN CLN12Did marked by the arrow; the right panel shows 
the difference image after template subtraction. 

(A color version of this figure is available in the online journal.) 



Figure 4. Top: HST composite false-color image of SN CLA1 ITib from CLASH 
data. The red channel is an average of WFC3/IR F105W, F110W, F125W, 
F140W, and F160W data. The green channel combines ACS F606W, F625W, 
F775W, F814W, and F850LP observations, and the blue channel is from ACS 
F435W and F475W data. The location of the SN is marked with a white circle. 
Bottom: the left panel shows a 5" square cutout of the ACS F850LP image taken 
on UT 2011 January 4, with SN CLAllTib indicated by the arrow; the right 
panel shows the difference image after template subtraction. 

(A color version of this figure is available in the online journal.) 


galaxy, so we expect negligible host-galaxy contamination in 
the SN photometry. 

Ground-based Sloan Digital Sky Survey (SDSS) DR7 images 
(Abazajian et al. 2009) of the CLN12Did host, along with the 
HST data, show the CLN12Did host to be an early-type, likely 
elliptical, galaxy. The CLASH data yield photometry for the 
host galaxy of F606W = 22.69 =b 0.03 AB mag and F160W = 
19.434 ± 0.003 AB mag. We obtained a ~5 hr spectrum 
(6 x 2700 s and 1 x 2500 s) of the host galaxy using the Multi- 
Object Double Spectrographs (MODS; Pogge et al. 2010) on 
the Large Binocular Telescope (LBT), covering 6300-9000 A 
with 15-40 A resolution. We also used Gemini North GMOS 
to obtain a 4.5 hr (9 x 1800 s) spectrum of the galaxy, covering 
5220-9420 A, with a resolution of ~6 A. A composite spectrum 
from the Gemini and LBT observations is shown in the middle 
panel of Figure 2. We derived respective redshifts of 0.851 ± 
0.001 and 0.852 =b 0.002 from these data, via cross-correlation 
with absorption-line templates. This redshift is consistent with 
the photo-z derived for the galaxy (z p hot = 0.92+°q q 5 7 ). We adopt 
z = 0.851 ± 0.001 for SN CLN12Did. At this redshift, the 
4" separation from SN to host-galaxy center corresponds to 
a projected physical separation of ~30 kpc (modulo lensing 
magnification). 


23. SN CLAllTib 

SN CLAllTib was discovered in CLASH ACS and WFC3/ 
IR images of A3 83 (Abell et al. 1989) in 201 1 January. Figure 4 
shows the SN in a color image (located at a = 02 h 48 m 01 s .27, 
5 = — 03°33 / 16 / . / 9, J2000). Aperture photometry for the ACS 
observations of SN CLAllTib is presented in Table 3. The 
WFC3/IR observations of the SN were taken as part of the HST 
follow-up program GO-12360 (PI: S. Perlmutter). The ACS 
photometry is from difference images with host-galaxy subtrac- 
tion, but no template images are available for the WFC3/IR 
observations (F105W, F125W, and F160W), so the reported IR 
photometry is based on the direct images. As the SN is located 
in a spiral arm (see Figure 4), the host galaxy could contaminate 
aperture photometry for these observations. For this reason we 
report point-spread function (PSF) fitting photometry of the SN 
in the WFC3/IR observations from Dolphot (a modified ver- 
sion of HSTphot; Dolphin 2000). Photometry of the host galaxy 
gives F606W = 22.88 =b 0.01 AB mag and F160W = 21.51 ± 
0.01 AB mag. 

We obtained Southern African Large Telescope (SALT) + 
Robert Stobie Spectrograph (RSS; Nordsieck 2012; Crawford 
et al. 2010) spectroscopy of the host of SN CLAllTib with a 
1200 s exposure using the PG0900 grating (range 6300-9300 A, 
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Table 3 

SN CLAllTib Photometry 


UT Date 

MJD 

Instrument 

Filter 

Exposure Time 

(s) 

Magnitude 

2010 Nov 18 

55518.91 

ACS 

F625W 

1032.0 

>27.206 

2010 Nov 18 

55518.92 

ACS 

F850LP 

1014.0 

>25.768 

2010 Nov 18 

55518.99 

ACS 

F775W 

1010.0 

>26.792 

2010 Dec 8 

55538.41 

ACS 

F850LP 

1032.0 

>25.802 

2010 Dec 8 

55538.43 

ACS 

F814W 

1060.0 

>26.872 

2010 Dec 28 

55558.45 

ACS 

F435W 

1032.0 

>27.385 

2010 Dec 28 

55558.46 

ACS 

F814W 

1092.0 

24.699 ± 0.049 

2011 Jan 4 

55565.97 

ACS 

F625W 

1032.0 

25.349 ± 0.074 

2011 Jan 4 

55565.98 

ACS 

F850LP 

1092.0 

24.251 ±0.091 

2011 Jan 18 

55579.35 

ACS 

F606W 

1032.0 

26.364 ±0.104 

2011 Jan 18 

55579.36 

ACS 

F814W 

1059.0 

24.394 ± 0.038 

2011 Jan 22 

55583.41 

ACS 

F775W 

1032.0 

24.762 ± 0.066 

2011 Jan 22 

55583.43 

ACS 

F606W 

1073.0 

26.551 ±0.127 

2011 Jan 24 

55585.08 

WFC3/IR 

F105W 

806.0 

23.755 ± 0.037 

2011 Jan 24 

55585.09 

WFC3/IR 

F125W 

806.0 

23.402 ± 0.036 

2011 Jan 24 

55585.14 

WFC3/IR 

F160W 

906.0 

23.357 ± 0.059 

2011 Feb 7 

55599.38 

ACS 

F814W 

1032.0 

25.082 ± 0.076 

2011 Feb 7 

55599.40 

ACS 

F435W 

1093.0 

>27.336 

2011 Feb 21 

55613.16 

ACS 

F850LP 

1994.0 

25.247 ±0.188 

2011 Mar 1 

55621.43 

ACS 

F850LP 

1076.0 

>25.677 


Notes. All SN photometry is reported in Vega magnitudes. The nondetections are listed as 3 a upper limits. Dolphot was used to derive PSF photometry 
for the WFC3/IR observations (see text). 


resolution 6 A), and detect a strong emission line (but close to a 
night-sky line) that we identify as [O n] at z = 1.144 =b 0.001. 
We confirmed this identification with the Visible MultiObject 
Spectrograph (VIMOS; Le Fevre et al. 2003) on the Very Large 
Telescope (VLT) observations in a ~4 hr spectrum of the host 
galaxy (range 3900-6700 A with 28 A resolution). The spectrum 
shows Mg ii and Fe n absorption, and [C m] A, 1909 emission, 
all at a consistent redshift of 1.143 ±0.001, which we adopt 
for SN CLAllTib. This redshift is consistent with the photo-z 
derived for the galaxy (z p hot = 1.05 ± 0.10). At this redshift, 
SN CLAllTib is projected ~25 kpc from the center of its host 
galaxy (uncorrected for lensing magnification). 

3. SUPERNOVA CLASSIFICATION AND 
LIGHT-CURVE FITS 

Classification of high-redshift (z 1) SNe is challenging 
(e.g., Rodney et al. 2012; Rubin et al. 2013; Jones et al. 2013). 
For example, a typical SN la has a peak luminosity in the rest- 
frame B band of M B & —19.4 mag. At z = 1.5 this becomes a 
T-band apparent magnitude of ~25.8. Key spectral absorption 
features at rest-frame optical wavelengths (such as the Si ii 
feature observed in SN la at ~6150 A) are redshifted into the IR, 
making them inaccessible to most ground-based spectrographs. 
With these constraints, definitive spectroscopic classification of 
high-z SNe requires prohibitively large investments of HST time 
(Jones et al. 2013). 

For this reason, we rely on photometric classifications for 
these three lensed SNe. We utilize a Bayesian photometric 
SN classifier called STARDUST: Supernova Taxonomy And 
Redshift Determination Using SuperNova ANAlysis (SNANA) 
Templates. An early version of this classifier was first intro- 
duced by Jones et al. (2013), and the complete version used 
for this work is presented in more detail by Graur et al. (2014) 
and S. A. Rodney et al. (in preparation). In brief, observations 
are compared against simulated SNe la and CC SNe, gener- 


ated using version 10.27k of the SNANA software package 36 
(Kessler et al. 2009b). For SNe la the STARDUST classifier 
uses the SALT2 light-curve model (Guy et al. 2010) extended to 
cover the rest-frame near-UV and IR wavelengths (see below), 
whereas for CC SNe, the classifier draws on the SNANA library 
of 16 SN Ib/c and 26 SN II templates, to be described in further 
detail by Sako et al. (2014). An important capability of STAR- 
DUST that we employ for this work is to allow an achromatic 
flux scaling factor for every model realization. This means that 
our STARDUST classifications rely only on the observed SN 
colors and light-curve shape, and do not depend on the intrinsic 
SN luminosities or our choice of cosmological model. This is 
especially important for fitting lensed SNe, as it allows us to 
classify the SNe independent of the lensing magnification. 

For additional verification of our photometric SN classifica- 
tions, we turned to the Photometric SuperNova IDentification 
software (PSNID; Sako et al. 2008, 2011) to classify the three 
SNe. PSNID is similar to the STARDUST classifier, as both use 
SALT2 to simulate SNe la and rely on the same 16 SN Ib/c and 
26 SN II templates to simulate CC SNe. However, PSNID has 
the advantage that it has been thoroughly tested, obtaining the 
highest figure of merit in the SN photometric classification chal- 
lenge (Kessler et al. 2010), and it is already publicly available as 
a component of SNANA. A downside of PSNID for this work 
is that it uses absolute magnitude priors in its classifications, 
which is not ideal for analyzing lensed SNe. 

For SNe la, we also employed a modified version of the Guy 
et al. (2010) SALT2 model to fit light curves and derive stan- 
dardized luminosity distances. The version of SALT2 presented 
by Guy et al. (2010) has been robustly tested, and was used 
in SNLS cosmology papers (Sullivan et al. 2011; Conley et al. 
201 1). To accommodate the wide wavelength range of our data, 
we extended the S ALT2 model to cover rest-frame near-UV and 
IR wavelengths (by default, the SALT2 fitter can only fit rest 


36 http://sdssdp62.fnal.gov/sdsssn/SNANA-PUBLIC/ 
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wavelengths of 2800-7000 A). For all three SNe, this extended 
SALT2 model produces similar results (within la) to those of 
the original Guy et al. (2010) model. In deriving distances, we 
converted our S ALT2 parameters representing light-curve shape 
and color excess (x 1 and c) to their respective SiFTO (Conley 
et al. 2008) parameters (s and C), using the conversion equa- 
tions from Guy et al. (2010). This was done in order to adopt 
the constants M, a, and p, from Sullivan et al. (2011) and also 
used by Jones et al. (2013). Specifically, the distance modulus 
was calculated using the following equation: 37 

dm SA LT2 =m* B ~ M + ot(s - 1) - pc, (1) 

with M = -19.12 =b 0.03, a = 1.367 zb 0.086, and p = 
3.179 zb 0.101. To match the (for our purposes, arbitrary) 
normalization of the SALT2 fitter used by Guy et al. (2010) 
to our results from SNANA, we applied an offset of 0.27 mag 
to the value of m* B reported by SNANA. 

We also used the Jha et al. (2007) version of MLCS2k2 
included in SNANA to fit the SNe la. To check for consistency, 
we compared the best-fit light-curve shape parameter A from 
MLCS2k2 to v 1 from S ALT2 using the relationship in Appendix 
G given by Kessler et al. (2009a). To compare the SALT2 color 
parameter c with MLCS2k2 Ay, we employ a linear relationship 
derived from Figure 42 of Kessler et al. (2009a), specifically 

c = (0.464 zb 0.021)Ay - (0.121 zb 0.014). (2) 

Finally, to ensure a consistent zeropoint between the SALT2 
and MLCS2k2 distance moduli, we applied both methods to 
a sample of SDSS Supernova Survey SNe la from Holtzman 
et al. (2008) and Kessler et al. (2009a). We added a zeropoint 
correction of 0.20 mag to the distance moduli reported by 
SNANA in the MLCS2k2 fit (dm M LCS2k2), to put them on the 
same scale as the SALT2 distances (dmsAm)- 

Although we report magnitudes in Tables 1-3, both the 
classification and SN la light-curve fitting were done in flux 
space. The upper limits for each SN were included in the fits 
with the measured flux and la uncertainties. 

3.1. SNCLOUCar 

SN CL012Car was located in a red disky galaxy, in which 
both SNe la and CC SNe could be expected. The STARDUST 
classifier prefers an SN la fit with probability (Pi a ) of 0.91. 
The probabilities for SN Ib/c (Pjb/c) and SN II (P n ) models 
are 10 -8 and 0.09, respectively. The best-fit SN la, SN Ib/c, 
and SN II models had x 2 / v (v = 11) of 1.12, 4.56, and 1.76, 
respectively. PSNID also favors an SN la fit at a high confidence, 
with Fi a = 0.99 and abest-fit x 2 /v = 1.06 (v = ll) 38 . We thus 
conclude that SN CL012Car was an SN la. 

Both MLCS2k2 and S ALT2 produce good light-curve fits for 
SN CL012Car. We show the SALT2 light-curve fit along with 
the la model errors in Figure 5. The SALT2 fit had a x 2 / v — 
13.3/11 = 1.21, and the fit parameters (xl = 1.375 zb 0.829 and 
c = 0.254 zb 0.049) are in the range of normal SNe la (Guy et al. 
2010). From these parameters, we derive a distance modulus of 
dm S ALT2 = 43.83 zb 0.24 mag for the SN. 


37 We use the notation “dm” for distance modulus, rather than the conventional 
/ji, to avoid confusion with the lensing magnification /jl. As usual, it is defined 
as dm = 5 log 10 di +25, where di is the luminosity distance in Mpc. 

38 PSNID uses a luminosity prior, as discussed above. Correcting the light 
curve for lensing, assuming reasonable values of magnification between 1 mag 
of demagnification and 3 mag of magnification, changes the Pi a value slightly 
(ranging from 0.87 to 0.99). 



Figure 5. SALT2 light-curve fit of SN CL012Car. The light-curve parameters 
xl and c are typical of normal SNe la, and the data are well matched by the 
model, with x 2 /v = 13.3/11. 

(A color version of this figure is available in the online journal.) 

The MLCS2k2 light-curve fit had a slightly worse, but still 
acceptable, goodness-of-fit (x 2 /v = 18.1/11 = 1.64). The 
model parameters were A = 0.007 zb 0.314 and A v = 0.635 zb 
0. 190 mag. The light-curve shape parameters from the two fitters 
agree at <la; specifically, the SALT2 xl value corresponds to 
A = —0.380 zb 0.160. The best-fit MLCS2k2 A v corresponds to 
c = 0.174 =b 0.090, which is also <la from the best-fit SALT2 
color. Moreover, we find that dniMLCS 2 k 2 = 43.71 zb0.16 mag 
(after applying the zeropoint correction), again consistent with 
the inferred dm S ALT2 at <la. 

For comparison, at z = 1.281, the distance modulus expected 
from ACDM with our adopted cosmological parameters is 
dmACDM = 44.76 mag, significantly higher than what is 
measured from the SN light curve. We discuss the implications 
of this finding in Section 4. 

3.2. SN CLN12Did 

SN CLN12Did was located in an elliptical galaxy, already 
indicating it was most likely an SN la (Cappellaro et al. 1999; 
Mannucci 2005; Foley & Mandel 2013). The STARDUST 
classifier prefers an SN la fit at high confidence (>5a). PSNID 
also favors an SN la fit with high significance (Pi a > 0.9999), 
with best-fit x 2 / v — 0-5 (y = 20). These consistent results 
yield the strong conclusion that SN CLN12Did was an SN la. 

We show the SALT2 SN la light-curve fit for SN CLN12Did 
in Figure 6. Similar to Figure 5, we display the light-curve points 
and the model curves with la errors. Upper limits on F435W 
and F475W nondetections are excluded from the plot for clarity, 
but all of the observations listed in Table 2 are used in the light- 
curve fit. The SALT2 parameters fall within the normal SN la 
range, with xl = 0.383 zb 0.208, c = —0.042 zb 0.021, and 
X 2 /v = 28.5/20 = 1.42. From the light-curve parameters, we 
derive dmsALT 2 = 43.33 zb 0.15 mag. 

MLCS2k2 produces a slightly better light-curve fit, with 
X 2 /v = 23.6/20 = 1.18. The best-fit parameters are A v = 
0.042 zb 0.059 mag and A = 0.055 zb 0.098. Using the 
conversion equations, the MLCS2k2 fit A v corresponds to 
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Figure 6. SALT2 light-curve fit of SN CLN12Did. Upper limits on F435W and F475W nondetections are excluded from the plot for clarity, but all of the observations 
listed in Table 2 are used in the light-curve fit. This model is a reasonably good match to the data, yielding x 2 / v = 1-4. The best-fit light-curve parameters are typical 
for a normal SN la. 

(A color version of this figure is available in the online journal.) 


SALT2 c = —0.102 d= 0.031, and the SALT2 xl predicts 
MLCS2k2 A = —0.212 zb 0.040. This is good agreement in the 
color parameter, but a slight (~2.5 cr) discrepancy in light-curve 
shape. The inferred distance is dm M LCS2k2 = 43.28 =b 0.09 mag, 
consistent with the SALT2 result. For comparison, at z = 0.85 1 , 
the cosmological parameters yield dm/\cDM = 43.65 mag. 

3.3. SN C LA 11 Tib 

SN CLN12Tib is located in a spiral arm of its blue, late-type 
host galaxy, an environment in which all SN types can be found 
(Cappellaro et al. 1999; Mannucci 2005; Li et al. 2011; Foley & 
Mandel 2013). The STARDUST classifier prefers a CC fit for 
this SN, returning probabilities Pu = 0.41 and P^/ c = 0.52, 
compared to the SN la fit probability F Ia = 0.07. 

We show the light-curve fit for SN CLAllTib in Figure 7, 
with each row giving the maximum-likelihood fit for a single 
SN subclass. The x 2 / v ( v — 14) values for the best-fit SN la, SN 
Ib/c, and SN II models are 3.81, 3.42, and 3.55, respectively. 
None of the three SN types are a good fit to the light curve; 
moreover, none of the templates provide a good fit to all the 
WFC3/IR data. There are only 42 CC SNe templates, and 
we do not expect them to represent the full range of SNe II 
and SNe Ibc. The SN II template that provides the best match 
to the observations is SDSS-018793 (SN 2007og), an SN IIP 
at z = 0.20 (Sako et al. 2014). For the SNe Ib/c, the best 
match comes from SN SDSS-017548 (SN 2007ms), an SN Ic at 
z = 0.0393 (Ostman et al. 2011). The sole HST WFC3/IR 
F125W observation of SN CLAllTib is close to the peak 
for both the SN IIP and SN Ic models; at z = 1.143, the 
F125W filter corresponds to approximately the rest V band. 
To further highlight the difficulty in photometrically classifying 
SN CLA1 ITib, in Figure 8 we show color-color diagrams from 
the ACS and WFC3/IR photometry taken during 2011 January 
22-24, where the SN straddles the typical colors of SNe la, Ib/c, 
and II, though favoring CC SN models. The full STARDUST 
classification uses these types of constraints with all available 
epochs. 


Our PSNID results disagree with those from the STARDUST 
classifier. PSNID favors an SN la fit for SN CLAllTib at 
greater than 4cr. However, the best-fit PSNID SN la model 
does not match the light curve well ( x 2 / v — 2.83, with v = 
14). To check these results, we fit SN CLAllTib with SALT2 
and MLCS2k2 directly using SNANA. The best SALT2 fit has 
normal light-curve parameters, with xl = 0.114 d= 0.486 and 
c = 0.048 zb 0.037, yielding dm S ALT2 = 44.07 d= 0.19 mag (for 
comparison, dm/\cDM = 44.45 mag). However, as with PSNID, 
the SNANA SALT2 fit requires an unacceptably large x 2 / v = 
57.8/14 = 4.13. These results do not change substantially 
if we use the restricted- wavelength version of SALT2; the 
fit remains poor. The MLCS2k2 results were much better 
(X 2 /v = 29.8/14 = 2.13), with best-fit light-curve parameters 
of = 0.389 =b 0.193 mag, A = -0.011 zb 0.244, and 
dm M LCS2k2 = 43.97zb0. 14 mag. Using the conversion equations, 
the MLCS2k2 fit Ay corresponds to SALT2 c = 0.059 zb 0.091, 
and the SALT2 xl predicts MLCS2k2 A = -0.201 zb 0.093. 
Interestingly, the MLCS2k2 and SALT2 parameters are similar, 
but the SALT2 fit is significantly worse. 

To summarize, our classification analysis of SN CLA1 ITib is 
inconclusive. The light curve was best fit by CC SN templates 
with STARDUST, while PSNID preferred an SN la fit. Similar 
to SN CL012Car, we corrected the light curve for a varying 
amount of lensing (from 1 mag of demagnification to 3 mag of 
magnification) and fit with PSNID to see if the luminosity prior 
was driving the classification. The value of Pi a varied greatly, 
with the probability dropping to 0 with magnifications >1.5 
mag. Because of these conflicting results, we further explore 
both CC SN and SN la models in our lensing analysis of SN 
CLAllTib (Section 4.3). 

4. GRAVITATIONAL LENSING MAGNIFICATION 

In this section, we derive the magnification of each SN in 
two independent ways. First, we use the SNe themselves, from 
the observed SN brightness and light-curve fits described in 
Section 3. Second, we independently determine the expected 
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Figure 7. STARDUST light-curve fits to CLAllTib for the three primary SN subclasses. The top row shows the maximum likelihood fit from an SN la (SALT2) 
model, the middle row displays the best-fit SN Ib/c model, and the bottom row shows the preferred SN II model. Optical data from ACS appear in the left column, 
while WFC3/IR observations are on the right. The CC SN models are labeled by the best-fit template SN, and all three models show the best-fit parameters. The 
downward pointing arrows indicate 3cr upper limits on the SN flux. The Am term designates the offset (mag) needed to make the template match the data, discussed 
further in Section 4.3. 

(A color version of this figure is available in the online journal.) 


magnification from cluster mass maps derived by the CLASH 
team from strong+weak-lensing features. 

One approach to inferring the SN magnification is simply 
measuring the offset between the inferred distance modulus and 
that predicted by the cosmological model (dm/\cDM)- In that 
case, however, statistical and systematic uncertainties in the 
cosmological parameters (which are themselves derived in part 
from SNe la over a wide redshift range) need to be included in the 
analysis. A more direct approach is to compare the CLASH SNe 
directly to a “field” sample at similar redshifts. For a sufficiently 
small redshift range, the relation between distance modulus and 
redshift can be taken to be linear. The field sample of SNe la 
define the zeropoint and slope of that relationship, and we derive 
the inferred magnification of our SN la relative to that locus, with 
all of these SNe analyzed in an identical way. This differential 
approach is insensitive to zeropoint offsets between the fitters 
and cosmological parameters, and only makes the assumption 
that all SNe la at nearly the same redshift are consistently 


standardizable (so it is less susceptible to evolution or other 
SN la systematics that vary with redshift). 

For the predicted magnification from the CLASH cluster 
mass maps, we analyzed both the strong- and weak-lensing 
features in the clusters. The identification of multiple images 
in the cluster fields was done using the 16-band HST imaging 
obtained in the context of the CLASH campaign. The Zitrin 
et al. (2009) method revealed nine multiple image systems 
in the case of A3 83 and seven systems in MACS J1720. 
Unfortunately, no multiple-image systems were found in the 
field of RX J1532, the cluster with the lowest total mass in 
this study. Where available, spectroscopic redshifts from the 
literature or the CLASH- VLT program (Balestra et al. 2013) 
determined the distance to the strong-lensing features. The 
reliable 16-band photometric redshifts (Jouvel et al. 2014) from 
the CLASH pipeline delivered the redshifts of systems without 
spectroscopic confirmation. In the case of weak lensing, we used 
5 -band CLASH Subaru imaging to derive weak-lensing shear 
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Figure 8. Color-color diagrams for SN CLAllTib photometry from MJD 55583-55585, illustrating the difficulty in photometric classification of this object. We 
utilize color-based classification so as to be unaffected by achromatic gravitational lensing. The data favor a CC origin for this SN (either SN Ib/c or SN II). The 
contours shown are derived from SNANA Monte Carlo simulations with 5000 simulated SNe for each of the three subclasses (so they do not reflect the relative 
frequency across classes). Shaded regions enclose 68% of each population and solid lines enclose 95%. 

(A color version of this figure is available in the online journal.) 


catalogs. The underlying pipeline for image reduction, shape 
measurement, and background selection is thoroughly outlined 
by Umetsu et al. (2012) and Medezinski et al. (2013). 

To combine the constraints from the weak and strong-lensing 
regime into a consistent lens model, we used SaWLens (Merten 
et al. 2009, 2011; Meneghetti et al. 2010; Umetsu et al. 
2012; Medezinski et al. 2013). This nonparametric method 
also reconstructs the mass distribution in the regime far away 
from the strong-lensing features, or where no strong-lensing 
features are available, and both cases apply in our study. By 
nonparametric, we mean that no a priori assumptions on the mass 
distribution of the lens were made. The method reconstructs 
the lensing potential, the rescaled and line-of-sight integrated 
counterpart of the Newtonian potential, on an adaptively refined 
mesh using both weak- and strong-lensing inputs. In the case 
of RX J1532, where no strong-lensing features were found, 
we use only the shear catalogs to constrain the lens mass 
distribution. From the lensing potential maps we directly derive 
the magnification (e.g., Bartelmann 2010) at the redshift of 
interest and read off its value at the position of each SN, as 
described individually below. In order to assign error bars to 
these values, we run 1500 realizations of each lens model by 


Table 4 

Comparing SN Magnifications with Lensing Predictions 


SN Name 

Lensing from Light-curve Fits 



Am /X 


Am /X 


(SALT2) 

(MLCS2k2) 

(Lensing Maps) 

SN CL012Car 

0.91 ± 0.25 

1.06 ±0.17 

0.83 ±0.16 

SN CLN12Did 

0.24 ±0.15 

0.15 ±0.09 

0.28 ± 0.08 

SN CLAllTib 



0.43 ±0.11 


Notes. Comparison of lensing predictions from light-curve fits and lensing 
maps for the three SNe. All of the A ra M table entries have units of magnitudes. 
For both SN CL012Car and SN CLN12Did, there is good agreement (within 
1 <t) between both lensing methods. The classification of SN CLAllTib was 
inconclusive, making it difficult to determine a magnification from the light 
curve (see Section 4.3). 


bootstrapping the weak-lensing input catalogs and randomly 
sampling through the allowed redshift error range of each 
strong-lensing feature (Merten et al. 2011). 

The independent magnification estimates, from both methods 
above, are summarized in Table 4. We consider each object 
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Figure 9. Hubble diagrams for SN CL012Car. The top panel shows SN la distances fit with SALT2 and the bottom panel has MLCS2k2 fits. The plots contain 18 
field SNe la from Riess et al. (2007, in black) and Suzuki et al. (2012, in green) as a comparison sample. Both light-curve fitters show that SN CL012Car (the large 
blue point) was significantly brighter than typical SNe la at similar redshifts. The solid black lines with the displayed equations give the best-fit linear model to the 
comparison sample (empirically approximating a cosmological fit over this small redshift range without any dependence on cosmological parameters)with the red 
dashed lines representing the lcr uncertainties. From these data, we predict that a typical unlensed SN la at z = 1.281 would have dmsALT 2 = 44.74 ± 0.07 mag and 
dniMLCS2k2 = 44.77 ± 0.05 mag. Compared to the inferred distance modulus of SN CL012Car, we find that it was magnified by Am /x> sALT2 = 0.91 ± 0.25 mag or 
Af?V,MLCS2k2 = 1.06 ± 0.17 mag. Both light-curve fitters give a consistent and significant magnification (greater than unity). 

(A color version of this figure is available in the online journal.) 
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individually below. Because we could not conclusively classify 
SN CLA1 ITib, we consider both the SN la and CC SN models 
for it (Section 4.3). 

4.1. SNCL012Car 

We showed in Sections 2.1 and 3.1 that SN CL012Car was 
an SN la at z = 1.281. In Figure 9, we display two Hubble 
diagrams for the SN. The top plot contains SN la distances fit 
with S ALT2 and the bottom plot with MLCS2k2. Both plots also 
show a comparison sample of 18 field SNe la from Riess et al. 
(2007) and Suzuki et al. (2012) in the range 1.14 z 1.41, 
fit in exactly the same manner as SN CL012Car. Consistently 
with either light-curve fitter, we find that SN CL012Car was 
significantly brighter than the typical SNe la at a similar 
redshift. The solid black line in both panels represents the 
best linear fit to the field SN la sample, with the fit slope 
and intercept as shown. 39 From this analysis, we can estimate 
the lensing magnification as determined by the SN itself. The 
linear fits to the field sample predict dmsAmCz = 1.281) = 
44.74 =b 0.07 mag and dm ML cs 2 k 2 fe = 1.281) = 44.77 zb 0.05 
mag. With these values, we derive that SN CL012Car was 
magnified by Am /X) sALT2 = 0.91 zb 0.25 mag (corresponding to 
a magnification /xsalt2 = 2.31 zb 0.54) in the SALT2 model, 


39 Three of the Suzuki et al. (2012) SNe la plotted in Figure 9 were also 
behind galaxy clusters and potentially lensed. However, the object with the 
highest estimated magnification still only has /z = 1.07, corresponding to 
A = 0.07 mag, well within the SN la scatter. Removing these three SNe la 
from the comparison sample has a negligible effect on our results. 


or Am lltjM LCS 2 k 2 = 1-06 ± 0.17 mag (/x M LCS 2 k 2 = 2.65 ± 0.42) 
with MLCS2k2. 

The cluster lensing magnification map (derived as described 
above) around MACS J 1720 for a source at the SN redshift (z = 
1 .281) is shown in Figure 10. From the map, we derive a lensing 
magnification for SN CL012Car of /x = 2.15 zb 0.33, which 
corresponds to A = 0.83 zb 0.16 mag. The magnification 
from the light-curve fits and the lensing maps agree to well 
within lcr. Moreover, SN CL012Car is inconsistent with being 
unlensed; its SN la fit requires A > 0 at ~4cr. 

4.2. SN CLN12Did 

We showed in Sections 2.2 and 3.2 that SN CLN12Did 
was an SN la at z = 0.851. In Figure 11, we display two 
Hubble diagrams for the SN, mirroring the analysis for SN 
CL012Car (see Section 4.1 and Figure 9). The comparison 
field sample contains 63 SNe la from the SNLS sample (Guy 
et al. 2010) in the range 0.75 z 0.95. Unlike the case for 
SN CL012Car, we see that SN CLN12Did falls well within the 
normal scatter for typical SNe la in the redshift range. Both 
light-curve fitters do find SN CLN12Did to be brighter than 
expected compared to the field sample, but not significantly so 
(~1 <t). Our inferred magnification estimates for SN CLN12Did 
are Am^sAm = 0.24 zb 0.15 mag (/x SA lt2 = 1.25 zb 0.18) and 
Aw /x, MLCS2k2 = 0.15 ± 0.09 mag (/XMLCS2k2 = 1.15 ± 0.10). 

The cluster lensing magnification map around RX J1532 for 
a source at z = 0.85 1 is shown in Figure 12. It predicts a lensing 
magnification for SN CLN12Did of A = 0.28 zb 0.08 mag 
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Figure 10. Magnification map for SN CL012Car (z = 1.281) behind MACSJ1720.2 + 3536 (z = 0.391). The image shows an RGB false-color background based 
on 12-band HST/CLASH optical and near-IR images of the cluster field. The location of the SN is marked with a white cross. The contours show the magnitude 
increase induced by the lensing magnification (A = 2.5 log 10 /x) of the cluster for a source at the SN redshift. The lensing magnification was derived from weak- 
and strong-lensing constraints jointly, and computed using the SaWLens lensing reconstruction algorithm (Merten et al. 2009, 2011). Multiple images used for the 
strong-lensing constraints in this system will be presented by A. Zitrin et al. (in preparation). 

(A color version of this figure is available in the online journal.) 



Figure 11. Hubble diagrams for SN CLN12Did. As in Figure 9, the top and bottom panels show SALT2 and MLCS2k2 fits, respectively. The field comparison sample 
consists of 63 SNe la from SNLS (Guy et al. 2010). Both light-curve fitters suggest slight (but insignificant) magnification for the SN, with Am^sALTC = 0.24 d= 0.15 
mag and Am M; MLCS 2 k 2 = 0.15 ± 0.09 mag. 

(A color version of this figure is available in the online journal.) 
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Figure 12. Magnification map for SN CLN12Did (z = 0.851) behind RXJ1532.9 + 3021 (z = 0.345). Similar to Figure 10, the color image is based on CLASH optical 
and near-IR images, with the SN location marked by a white cross and lensing magnification (A m^, in mag units) shown with red contours. Because RX J1532 is a 
much weaker lensing cluster than either MACS J1720 or A383, this map was derived using weak-lensing constraints alone, computed with the SaWLens algorithm 
(Merten et al. 2009, 2011). The resolution of the reconstructed map from weak-lensing constraints is much poorer than for the other two clusters, so the displayed 
contours are heavily smoothed, and the Monte Carlo magnification estimate at the SN position is actually A = 0.28 ± 0.08 mag. 

(A color version of this figure is available in the online journal.) 


(/x = 1.29 ± 0.09). Just like in the case of SN CL012Car, 
the lensing map prediction is consistent with the magnification 
inferred from the SN light curve, though for SN CLN12Did the 
magnification is only detected at low significance. 

4.3. SN C LA 11 Tib 

We were unable to classify SN CLAllTib conclusively (see 
Section 3.3), so we take a slightly different approach to the 
analysis for this object. Figure 13 shows the predicted lensing 
magnification map for a source at the redshift of SN CLA1 ITib 
(z = 1.143) behind A383. From the map, we derive a lensing 
magnification for the SN of /x = 1. 48 ±0.15, which corresponds 
to A/77 fl = 0.43 ±0.11 mag. 

If we assume that SN CLA1 ITib was an SN la, repeating the 
analysis above (as for SN CL012Car and SN CLN 12Did) yields 
magnification estimates for SN CLAllTib of Am^sAun = 
0.52 ± 0.20 mag and Am /XjM LCS2k2 = 0.64 ± 0.15 mag. 
However, as mentioned in Section 3.3, the SALT2 light-curve 
fit is unacceptable, with x 2 /v = 57.8/14 = 4.13, making 
the inferred distance unreliable. The MLCS2k2 light-curve fit 
(X 2 / v = 29.8/14 = 2.13) is much better, and agrees well 
with the lensing map prediction (within 1 a). However, due to 
the discrepancy of the goodness of fit between MLCS2k2 and 
SALT2, we hesitate to definitively classify the SN as an SN la. 

To compare both SN la and CC SN models in the lens- 
ing analysis for SN CLAllTib, we show in Figure 14 the 
goodness-of-fit for both sets of models. The plot gives the 


best-fit SN la results as above, as well as SN II and SN lb/ 
c templates (restricted to those with x 2 / v < 5). For the CC 
SN templates, STARDUST calculates a magnitude offset Am 
(shown also in Figure 7) that is required for the best match of 
the template to the data, and we show this on the abscissa of 
Figure 14. 

We found that the best-fitting CC SN models (albeit at 
X 2 /v > 3) from STARDUST were the SN Ic SDSS-017548 
and the SN IIP SDSS-018793. The SN Ic and SN IIP models 
had best-fit model m j = 23.52 ± 0.08 mag and m j = 
23.48 ± 0.08 mag, respectively. The observed red color of SN 
CLA1 ITib matches the SN Ic model with no relative extinction, 
but the SN II model requires A v = 0.7 mag. Including this, we 
find the unextinguished absolute magnitudes of the models are 
M r — — 19.34±0.08mag(SNIc)or Mr — —20.43 ±0.08 mag 
(SN IIP). These are quite bright for CC SNe, though some 
of this could be a result of the lensing magnification. The 
peak Mr required by both models is much brighter than the 
original templates; the SN SDSS-017548 and SN SDSS-018793 
spectral energy distributions have Mr = —18.00 mag and 
Mr = —18.55 mag, respectively. The offset between the 
template and data (Am; Figures 7 and 14) is —1.16 mag 
for the SN Ic model and —1.18 mag for the SN IIP model. 
These would be our best estimates for the lensing magnification 
based on the SN light curve, but an abundance of caution is 
warranted, as these estimates assume that SN CLAllTib had 
the same intrinsic luminosity as the best-fit template, a dubious 
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Figure 13. Magnification map for SN CLAllTib (z = 1.143) behind A3 83, similar to Figures 10 and 12, computed from the SaWLens algorithm (Merten et al. 2009, 
2011) from strong and weak-lensing constraints. The strong lens multiple images in this system were presented by Zitrin et al. (2011). 

(A color version of this figure is available in the online journal.) 
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Figure 14. Reduced y 2 (with 14 degrees of freedom) vs. magnitude offset for SN CLA1 ITib. The magnitude offset is a proxy for A from the light-curve fit model, 
as described in the text. The cumulative probability for a given xl value is displayed on the right-hand side. The blue shaded region corresponds to the lcr range of 
A derived from the lensing map in Figure 13. The black and blue points correspond to SN Ib/c and SN II models, respectively, while the red (SALT2) and magenta 
(MLCS2k2) points are the best-fit SN la models. Only the MLCS2k2 fit is consistent with both the SN light curve and the lensing map magnification prediction, but 
the CC SN templates do not span the full range of CC SN luminosities and light curves. 

(A color version of this figure is available in the online journal.) 
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proposition given the large luminosity scatter of CC SNe (Li 5. DISCUSSION AND CONCLUSION 

et al. 2011; Drout et al. 2011; Kiewe et al. 2012). Because of 
this, if we assume a CC SN fit, we cannot claim that there is any 
consistency or inconsistency between the lensing prediction for 
SN CLAllTib and its brightness. 
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SN CL012Car clearly represents our most exciting result. 
We find that it is likely an SN la at z = 1.281, gravitationally 
magnified by ~1.0 =b 0.2 mag. The magnification derived from 
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two independent methods (the SN light curve and the cluster 
lensing maps), were consistent to within 1 a . SN CL012Car 
is the first SN la which is both measurably magnified (i.e., 
with A > 0 at high significance) and for which the lensing 
magnification can be precisely and independently derived. SN 
CLN12Did is also an SN la with a consistent (but not significant) 
SN-based and model-based lensing magnification, whereas the 
unclear classification of SN CLA1 ITib does not allow for a test 
of the lens model prediction. 

The lens models used here were constructed jointly with both 
strong-lensing constraints (multiple images) and weak-lensing 
constraints. In principle, the strong-lensing multiple images 
directly constrain the deflection field, which is the gradient of the 
potential, while the weak-lensing ellipticity measurements fur- 
ther out constrain the (reduced) shear, a combination of second- 
order derivatives of the potential. As an alternative to the analy- 
sis presented here (where we independently tested the estimated 
magnifications), lensed standard candles like SNe la can play 
an important role as an additional direct , local constraint on the 
magnification (which is a function of the convergence k, and 
the shear y , both constructed from second-order derivatives of 
the potential) in the lens model (e.g., Riehm et al. 201 1). 

Gathering large samples of lensed SNe la is clearly of interest. 
Due to the magnification power of the foreground lenses, lensed 
SNe could be expected to be seen to higher redshifts than 
field SNe (for a given limiting magnitude), with applications 
to extending the redshift coverage of cosmological parameter 
measurements (if the lensing magnification can be estimated 
precisely and accurately enough, e.g., Zitrin et al. 2013) or SN 
rate measurements to constrain progenitor models (e.g., Graur 
et al. 2014; Rodney et al. 2014). There is also the exciting 
prospect of finding multiply imaged SNe, with measured time 
delays as a new cosmographic constraint (e.g., Holz 2001; 
Goobar et al. 2002; Oguri & Kawano 2003; Oguri et al. 2003). 

SN la luminosity distances, combined with lens model 
magnification estimates (which depend on angular diameter 
distances to and between the lens and source) or time-delay 
distances from multiple SN images, could also open up new tests 
of general relativity (e.g., Daniel & Linder 2010; Jain & Khoury 
2010; Schwab et al. 2010) or even fundamental cosmological 
assumptions, like the distance duality relation (e.g., Lampeitl 
et al. 2010; Liang et al. 2013). 

Though the CLASH survey has ended, the three lensed SNe 
found behind CLASH clusters herald the promise of many more 
such discoveries in the near future, including from wide-held 
ground-based surveys like the Dark Energy Survey (Frieman 
& Dark Energy Survey Collaboration 2013; Bernstein et al. 
2012), Pan-STARRS (Rest et al. 2013 and references therein), 
and ultimately LSST (LSST Science Collaboration et al. 2009); 
ground-based surveys that target massive galaxy clusters, such 
as CluLeSS (S. W. Jha et al., in preparation) or SIRCLS 
(Graham et al. 2013); and targeted space-based surveys like 
the HST Frontier Fields (PI: M. Mountain; with SN follow- 
up observations in program GO- 13396, PI: S. Rodney). SNe 
la were the key to the discovery of the accelerating universe; 
perhaps gravitationally lensed SNe la will play a starring role 
in further illuminating the dark universe, probing not just dark 
energy, but dark matter as well. 
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